arXiv:1504.00577vl [astro-ph.HE] 2 Apr 2015 


Submitted to ApJ March 7, 2015; accepted April 2, 2015 

Preprint typeset using style emulateapj v. 5/2/11 


RULING OUT IC/CMB X-RAYS IN PKS 0637-752 AND THE IMPLICATIONS FOR TEV EMISSION FROM 

LARGE-SCALE QUASAR JETS 

Eileen T. Meyer^ 

Space Telescope Science Institute, Baltimore, MD 21218, USA 

Markos Georganopoulos^ 

University of Maryland Baltimore County, Baltimore, MD 21250, USA 
William B. Sparks 

Space Telescope Science Institute, Baltimore, MD 21218, USA 


Leith Godfrey 

ASTRON, Netherlands Institute for Radio Astronomy, P.O. Box 2, 7990 AA Dwingeloo, Netherlands 


James E. J. Lovell 

School of Physical Sciences, University of Tasmania, Private Bag 37, Hobart, Tasmania 7001, Australia 

AND 

Eric Perlman 

Florida Institute of Technology, Melbourne, FL 32901, USA 
Submitted to ApJ March 7, 2015; accepted April 2, 2015 

ABSTRACT 

The Chandra X-ray observatory has discovered dozens of resolved, kiloparsec-scale jets associated 
with powerful quasars in which the X-ray fluxes are observed to be much higher than the expected level 
based on the radio-optical synchrotron spectrum. The most popular explanation for the anomalously 
high and hard X-ray fluxes is that these jets do not decelerate significantly by the kiloparsec scale, but 
rather remain highly relativistic (Lorentz factors P ^10). By adopting a small angle to the line-of- 
sight, the X-rays can thus be explained by inverse Compton upscattering of CMB photons (IC/CMB), 
where the observed emission is strongly Doppler boosted. Using over six years of Fermi monitoring 
data, we show that the expected hard, steady gamma-ray emission implied by the IC/CMB model 
is not seen in PKS 0637-752, the prototype jet for which this model was first proposed. IC/CMB 
emission is thus ruled o ut as the source of the X-rays, i oining recent results for the jets in 3C 273 
(using the same method; [Mever fc CeorganoDOulosl[20141 and PKS 1136-135 (using UV polarization; 
iCara et ^1201311 . We further show that the Fermi observations give an upper limit of J <6.5 for the 
four brightest X-ray knots of PKS 0637-752, and derive an updated limit of d <7.8 for knots A and 
B1 of 3C 273 (assuming equipartition). Finally, we discuss the fact that high levels of synchrotron 
X-ray emission in a slow jet will unavoidably lead to a level of angle-integrated TeV emission which 
exceeds that of the TeV BL Lac class. 

Subject headings: galaxies: jets — galaxies: active — quasars: individual (PKS 0637-752, 3C 273) 


1. INTRODUCTION 

In August 1999, the Chandra X-ray Observatory ob¬ 
served its first celestial target, quasar PK S 0637-752, dur¬ 
ing the initial focusing o f the telescope (iSchwartz et al.l 
120001 lOiart as et ^1200011 . Along with the bright quasar 
core, Chandra unexpectedly detected X-rays from the 
kiloparsec (kpc) scale relativistic jet (previously known 
from radio imaging, Figure 1). Unlike the synchrotron 
spectrum of lower-power FR 1 jets like M87 which eas¬ 
ily e xtend up to X-ray energies (e.g. iWilson fc Yand 
l2002fl . the synchrotron spectrum of powerful kpc-scale 
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jets (including PKS 0637-752) generally peak at or be¬ 
low the IR/Optical band. The X-rays detected in 
the kpc-scale jet of PKS 0637-752 were much brighter 
than expected from the radio-optical synchrotron spec¬ 
trum, or indeed from either synchrotron self-Compton 
(SSC) or inverse Compton upscattering of ambient CMB 
photons (IC/CMB) assuming mildly re lativistic kpc- 
scale flows under equipartition conditions (jChartas et al.l 
l2000fl . Further, the X-ray spectrum of the jet was re¬ 
markably hard, with a photon index of 1.76 ± O.I (Fig¬ 
ure [2]). 

Multi-epoch measurements of sub-parsec scale jets of 
powerful quasars with Very Long Baseline Interferome¬ 
try (VLBI) have detected superluminal proper motions 
which imply that these jets start out hi ghly relativis- 
tic, with Lorentz factors (P) of 10-50 (jJorstad et al.l 
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I2005t iLister et alJ I2009D . Though it had long been 
supposed based on population studies that jets decel¬ 
erate and are at most mil dly relativistic by the time 
they reach the kpc scale (e .g. lArshakian fc Longaiill2nn4 
iMullin fc Hardca stlell2009ll. no d i rect m eas urements have 
confir med this. iTavecchio et alJ (j200r)ll and iCelotti et alJ 
(|2001ll thus suggested that the X-rays from the jet in 
PKS 0637-752 could be explained by IC/CMB emission 
if the jet remained highly relativistic (T ^10), and was 
pointed at a fairly small angle to our line of sight (6°). 
This produces a much larger Doppler boosting factor 
[S '^10) and enables the IC/CMB X-ray flux to match 
the observations. 

Over the past decade and a half since the launch of 
Chandra, dozens more kpc-scale quasar jets with anoma- 


louslv hard and/or hi 
fe.g. [Sambnma et al. 

jh X-rav fluxes have been detected 
[20011 I 2 OO 2 I: [Siemiginowska et al. 

200,11: ISambruna et al.l 120041: IMarshall et al.l 

2005; 

Harris & Krawczvnskil [20061: [Siemiginowska et al.[ 

2007; 

Marshall et al.[ [20 1 It 

[Kharb et al.[ [20I2I: [Godfrey et al. 

20I2all. The IC/CMB model has been bv far the 


most popular explanation of thes e X-rays, t hough 
problems have been not ed (e.g iHardcastl^ I2006I 
[Harris fc Krawczvnskil [^00611 . Besides the unconfirmed 
fast speeds required on the kpc scale, IC/CMB often 
requires the jet to be pointed very close along our 
line-of-sight, leading to a deprojected jet length longer 
than 1 Mpc, comparable to, o r greater than the larges t 
known giant radio galaxies (jPermer fc Atovai] 1200411 . 
Further, the electrons responsible for upscattering the 
CMB into the Chandra band are at much lower energies 
than are traced by radio observations. This extension of 
the electron energy distribution is energetically costly, 
in some cases leading to ‘super-Eddi n gton’ jet power 
requi re ments (| Dermer fc Atova nI 12004 IJchivama et al.l 
1200611 . iCcorganoDOulos fc Kazana^ (120041 1 also noted 
that the commonly observed severe decrease in X-ray to 
radio flux ratio w ith distance down the jet (e.g. 3C 273, 
[Jester et al.ll2006ll can only be reconciled wit h IC/CMB 
X-ray s if the jet is strongly decelerating. [Hardens tI3 
(|2006ll . seeking to verify this for a sample of quasar jets, 
also noted that this explanation still leaves significant 
discrepancies with observations, namely the lack of 
similar deceleration profiles at kpc-scales for more 
misaligned jets. These problems lead to the suggestion 
that the X-rays could alternatively be synchrotron 
emission from a second elec tron population in the 
jet, albeit of unknown origin dA tovan fc Derme il[2004 
Harris et al.[[200l : [Kataoka fc Stawarz[[2005 jHardcastlel 


2006t [Jester et al.][2006l: [Uchivama et al.l[2006[). 

Despite a significant effort by the community to amass 
high-resolution radio, optical, and X-ray imaging of 
dozens of quasar jets, the fundamental problem up to 
now has been that fitting the radio-to-X-ray spectral 
energy distribution (SED) alone cannot distinguish be¬ 
tween t he IC/CMB and sy nchrotron explanations for the 
X-rays (|Cara et al.l[20I3h . The difference in power re¬ 
quirements between the two mechanisms is great, as is 
the extremely different idea of jet structure that they im¬ 
ply. Discriminating between these models is essential to 
make progress on the the physics of jets and their impact 
on their environments. 

It was with an eye to resolving t his lon g-standing dead¬ 
lock that [Georganopoulos et al.l ([200611 . hereafter G06, 
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Fig. 1.— Up 

per Panel: Chandra X-ray image of PKS 0637-752 

Chartas et al.l 

2000Q. with A'rCJA 18 GHz radio contours overlaid 

(iodfrev et alJ 

2Cn2al. Lower Panel: Optical image of PKS 0637- 


752, taken with ACS/WFC on HST in F475W HMehta et al.l[2?i(1^ 1 
with ATCA 18 GHz radio contours overlaid. 


suggested that Fermi Large Area Telescope (LAT) ob¬ 
servations could confirm or rule out the IC/CMB mecha¬ 
nism for the X-rays, by detecting (or not) the high level of 
gamma-ray emission this mechanism predicts (Figure [5]). 
We have previously looked for this gamma-ray emission 
from the jet of 3C 273, and ruled out IC/CMB gamma- 
rays from (the brightest) knot A alone at the >95% level, 
and f rom knots A through D1 combi ned at the >99.9% 
level (jMever fc Ceorganor)oulosll20l4 hereafter M14). 

In this paper we report new Fermi observations of 
PKS 0637-752 which show that the expected steady 
gamma-ray emission from the IC/CMB mechanism is 
also ruled out by deep upper limits at the 99.9% level. 
We also present updated limits for 3C 273, showing that 
the expected gamma-rays from IC/CMB are now ruled 
out at the 99.99% level in more than one Fermi energy 
band. We show that the deep upper limits at GeV ener¬ 
gies place interesting constraints on the Doppler beaming 
factors which implies that these jets are not highly rela¬ 
tivistic on the kpc scale. We further discuss the surpris¬ 
ing implications of slow, synchrotron X-ray jets on our 
understanding of the total radiative output of quasars, 
especially at TeV energies. 

2. METHODS 

2.1. The Fermi Test of IC/CMB 

As first noted by G06, the shape of the IC/GMB spec¬ 
trum is constrained to match the synchrotron spectrum, 
with a shift in frequency and luminosity solely deter¬ 
mined by the factor B/6 where S is the Doppler beaming 
factor and B the magnetic field strength. From G06, we 
have: 


__ 

eBS^"^/[2TTmec{l + z)] 
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Fig. 2.— A depiction of the Fermi test for the X-ray emission 
mechanism in large-scale jets. Here we show the radio, optical, and 
X-ray photometry for the knots (total) in the jet of PKS 0637-752 
as black points. The thin black line shows the synchrotron fit to 
the data. This synchrotron spectrum must be shifted by a factor 
3.4x10^® in frequency and 3.05x10^ in luminosity, corresponding 
to B/5 = 1.5x10“^ G, in order for the IC/CMB spectrum (thick 
black line) to match the X-ray flux level. As shown, this implies 
a high level of gamma-ray emission which should be detectable 
with Fermi. Alternatively, the X-rays could be due to a second 
synchrotron component (thin orange line). 


L, 327rUcMB{l + z)*S* 

Ls ^BS)^ ’ ^ ’ 

where Vc and Vg {Lc, Lg) are the observed Compton and 
synchrotron frequencies (luminosities) emitted by elec¬ 
trons of Lorentz factor 7 , e and me are the electron 
charge and mass, and zzcmb = 1.6x10^^ Hz is the CMB 
peak frequency at z = 0. However, if the observed X-ray 
fluxes are to be produced by the IC/CMB mechanism, 
then the value of B/5 is already uniquely determined by 
the requirement to match the X-ray flux level, at which 
point there is no freedom at all in the rest of the spec¬ 
trum. The peak of the IC/CMB spectrum will fall in the 
GeV band, as shown in Figure [5J Note that this predic¬ 
tion is not predicated on any particular (e.g., equiparti- 
tion) magnetic field strength. 

The Eer 77 i*/LAT lacks the spatial resolution to detect 
the jet separately from the gamma-ray bright quasar 
core, which is only 10" away - the Eermz/LAT 68 % 
containment radius is on the order of tenths of a de¬ 
gree to degrees. However, in powerful quasars the in¬ 
verse Compton core emission generally peaks at a few 
MeV, producing a soft and extremely variable spectrum 
in the Fermi band, with long periods of relative qui¬ 
escence. Indeed, PKS 0637- 752 was detected i n the 2nd 
Fermi source catalog ('2FGL. lNolan et aI1l2012ll as source 
2FGL J0635.5-7516, with a very soft spectrum (photon 
index of rp=2.71) and high degree of variability (vari¬ 
ability index = 347). In contrast, the IC/CMB emission 
from the large scale jet is expected to be harder and 
completely non-variable. The latter property allows us 
to combine the Fermi data taken only when the quasar 
core is in a low state to try to detect or place limits on 
the IC/CMB emission. 


2 . 2 . Fermi Analysis of PKS 0637-752 

We first combined the all-sky weekly LAT event and 
spacecraft files for weeks 9 through 325 of the Fermi 
mission, corresponding to Fermi Mission Elapsed Time 
(MET) from 239557417 to 430608212 and calendar dates 
2008 August 4 to 2014 August 24. In order to analyze the 
region around PKS 0637-752, we used the publicly avail¬ 
able ‘quickAnalysis’ scrip10 to run the Fermi analysis 
tools and generate the filtered event file, livetime cube, 
and exposure map, using a region of interest (ROI) of 
lO*^ and an otherwise default configuration. The start¬ 
ing source list was generated from the publicly available 
make2FGLxml script, which generates the xml file pre¬ 
populated with 2FGL catalog sources. We used a binned 
maximum likelihood to get an initial fit for all the catalog 
sources in our ROI. We also included sources a further 
5“^ out from our ROI, but with spectral parameters fixed 
to the 2FGL catalog values. PKS 0637-752 was detected 
with a very high test-statistic (TS, roughly significance 
squared) value of 289, a 100 MeV to 100 GeV photon 
flux of 3.16x10“® s“^ cm“^, and a photon index Pp = 
2.64, similar to the value reported in the 2FGL catalog. 

We checked for additional significant sources within 7° 
of PKS 0637-752 but not in the 2-year LAT catalog by 
making a TS residual map. The TS map was created by 
first freezing all source parameters in the initial model 
to their best-fit values, and then checking what TS value 
a test point source at different locations would be given. 
The resulting initial TS map is shown at left in Figure [3l 
In two areas the coarse 0.5*^ binning did not allow separa¬ 
tion of what appeared to be multiple components, and so 
a finer 0.125'^ binning was used over these smaller areas 
to determine the rough location of new sources. 

Rough starting positions of apparent new significant 
sources were measured from the TS map by hand, only 
considering as candidates those with a central pixel value 
(TS) > 20. Each potential new source was then localized 
one at a time as follows. A power-law point source was 
added at the rough location of the excess TS, with nor- 
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Fig. 3.— Left: An initial TS map of the region around PKS 0637- 
752, showing the excess TS present in 0.5° pixels over the best- 
fit likelihood model using 2FGL catalog sources. The large circle 
marks the position of PKS 0637—752. The smaller circles mark 
the positions of the six new sources not present in the 2FGL cata¬ 
log (corresponding to regions of excess TS with pixel values >20). 
Right: The updated TS map (same field-of-view and binning) 

after the six new sources were localized and fit with a binned like¬ 
lihood. 


^ The public scripts mentioned in the text are available at 
http: //fermi.gsfc.nasa.gov/ssc/data/analysis/user/ 
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TABLE 1 

New sources in ROI of Targets 


near PKS 0637-752 near 3C 273 


R.Aj2000 

(deg) 

Dec j 2000 
(deg) 

TS 

RAj2000 

(deg) 

Decj 2 ooo 

(deg) 

TS 

82.43046 

-72.74572 

57.5 

192.82676 

-2.00263 

70.7 

81.13040 

-69.60575 

60.5 

190.96827 

-2.29561 

82.3 

78.98804 

-72.72570 

28.6 

187.15935 

-3.29476 

60.2 

86.33797 

-70.34640 

39.6 

184.47577 

-0.48239 

92.5 

119.70430 

-80.70430 

9.9 

193.43690 

3.47391 

224.6 

118.83300 

-80.32970 

16.1 

192.63313 

2.25116 

65.7 


malization and powerlaw index free. Sources within 5° 
of the new source also had normalization/index parame¬ 
ters free. We ran a binned likelihood to get a starting fit 
for the new source plus other sources in the area. The 
ht parameters were then frozen again, so that we could 
next optimize the RA/Dec location of the new source. 
The existing tool, gtfindsrc, only works for unbinned 
likelihood analysis, so we built our own binned version of 
the tool which works in the same way. Using the frozen 
model, we used the python minimize function in the 
scipy package (L-BFGS-B method) to optimize the log- 
likelihood value versus the RA and Dec position, given 
a reasonable range of about 1 degree around the start¬ 
ing positions noted by hand. We then update our new 
source in the model with this optimized position (posi¬ 
tions always remain fixed when optimizing over spectral 
parameters). We un-freeze the normalization and index 
of the new source and the local sources and run a final 
binned likelihood optimization to get an updated model 
of all sources in the ROI. This is then the starting model 
for the next new source to be added, until all new sources 
have been localized and a final optimized fit derived. For 
the 7° ROI around PKS 0637-752, six new sources were 
added to the model, and an updated TS map run from 
this larger source list shows that the excess TS previously 
seen is now gone (right panel of Figure [3]). A list of the 
new sources with their location and TS value is given at 
left in Table [H 

Modified Julian Day 
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Fermi Mission Eiapsed Time (MET) [sec /10® ] 


Fig. 4. — Upper: Lightcurve of PKS 0637-752. The total 100 
MeV - 100 GeV photon flux for PKS 0637-752 in 10.5-day (total 
GTI) bins versus the mean Mission Elapsed Time (MET) of the 
bin. Upper limits are shown where TS<10. The gray vertical lines 
correspond to the start of a calendar year, beginning with 2009. 
Lower: The TS value of the bin versus MET. 


As the PKS 0637-752 core is a significant Fermi source, 
our approach to detecting and/or setting limits on the 
IC/CMB gamma-ray emission exploits the variability of 
the core which cannot be spatially resolved separately 
from the large-scale jet due to the poor angular resolution 
of Fermi. During times when the blazar is quiescent, the 
hard, steady emission from IC/CMB will either appear as 
a steady plateau, or else the upper limits generated will 
place constraints on the level of the IC/CMB emission. 

In order to build a lightcurve of the core, the full 6-year 
dataset was divided into bins of equal good time interval 
(CTI) time, totaling 10.5 days. We then used our up¬ 
dated (2FCL -I- 6 new) model described above and ran 
a binned likelihood to fit PKS 0637-752 as a power-law 
source, with sources more than 5° away hxed. The re¬ 
sulting lightcurve for the core over the full time range is 
shown in Figure SI with the 100 MeV - 100 CeV pho¬ 
ton flux shown on top and the corresponding TS shown 
below. The clear variability in the light curve indicates 
that the total Fermi flux is dominated by the core. 

We next began a ‘progressive binning’ analysis, in 
which the lightcurve bins were ordered from lowest to 
highest TS value. Of the entire set of 75 time bins, 27 
showed a TS level consistent with zero for the location 
of PKS 0637-752 (upper limits in lower panel of Fig¬ 
ure SD- Starting from these 27 bins combined, we pro¬ 
gressively combine the event files for the lowest bins plus 
the next lowest bin in TS, at each step optimizing the fit 
of PKS 0637-752 and the sources within 5° with a binned 
likelihood. We repeated adding the next-highest bin and 
getting the maximum likelihood fit until all bins had 
been added together. Note that this re-combining of the 
lightcurve in a discontinuous way is appropriate for de¬ 
riving a limit on the large-scale jet because the IC/CMB 
emission is predicted to be completely non-variable, and 
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Fig. 5.— The results of the progressive binning analysis on 
PKS 0637-752. The upper axis gives the flux while lower 
axis gives the total number of bins combined (where ordering is 
based on TS value and not time order), starting from the 27 bins 
with TS<0.01. The points give the upper limits (colored dots) or 
detected fluxes (TS>10, black triangles) for each of the five Fermi 
energy bands (red, orange, green, cyan, navy from lowest to high¬ 
est energy). Errors on the detected fluxes are only shown for the 
lowest-energy band for clarity but are similar across bands. For 
the highest-energy band (10 - 100 GeV) no detection is ever made. 
The increasing detected fluxes indicate that the quasar core is be¬ 
ing detected in the other bands. 
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Fig. 6.— Left: The SED for the large-scale jet of PKS 0637-752. Data for the four brightest X-ray detected knots combined is shown as 
black points. The SED of the X-ray brightest knot, wk8.9, is plotted with blue points. Right: The SED for the knots of 3C 273, with black 
points for knots A and B1 combined, and blue for X-ray brightest knot A only. For both plots, the Fermi 95% upper limits are shown in 
red. 


thus there is no risk of any selection effect via variabil¬ 
ity. The variable core clearly dominates the flux levels 
determining the ordering, and is disconnected from the 
jet in any case. At each step we evaluated the TS and 
flux level in the five canonical Fermi energy bands of 
0.1-0.3, 0.3-1, 1-3, 3-10 and 10-100 GeV, calculating the 
95% upper limit flux value when TS<10 in any given 
bin. Previous work on 3C 273 has shown that the exact 
method of ordering the bins (whether by using the TS 
value or the total flux or upper limit value for the bin) 
does not significantly affect the order or resulting upper 
limits (M14). 

As shown in Figure [HI the highest energy bands gave 
upper limits which decreased with the increasing expo¬ 
sure as more bins are added, where we have color-coded 
the upper limits in the 5 energy bands, and black trian¬ 
gles indicate a significant detection in the band. Note 
that in the lower-energy bins, where the PKS 0637-752 
quasar core dominates due to its soft spectrum, the upper 
limits reach a minimum rather quickly, and generally in¬ 
crease before becoming detections. It must be noted that 
the detected gamma-ray emission in these bands is from 
the quasar core, not the large-scale jet, based on the soft 
spectrum, and the fact that the emission level rises as 
more bins are added (showing that the source is indeed 
variable and that the bins are ordered by flux level). If 
we were to detect the large-scale jet, we would expect 
it to contribute significantly in the two highest-energy 
Fermi bands, and also to remain steady in flux level af¬ 
ter the exposure became long enough to detect it. This 
latter point is important for distinguishing a detection 
of IC/CMB gamma-rays versus the tail of the Compton 
emission from the core, which should (like the low-energy 
bins) rise in flux as more bins are added. While the fourth 
energy band detected points do not rise as quickly as the 
first three, the flux level is far above the upper limit de¬ 
rived after 50 bins, so cannot be the steady emission of 
the large-scale jet, which must be below this limit. The 
highest bin never shows a significant detection of either 
component. 


2.3. Fermi Analysis of 3C 273 


We re-analyzed the Fermi data for 3C 273 using the 
6-year dataset to compare with the results from M14 
using 4.5 years of data, as the core remained relatively 
quiescent over the additional time elapsed. We followed 
the same procedure as outlined above for PKS 0637- 
752, finding six new sources within 7° of the position of 
3C 273, listed at right in Table[T] The core of 3C 273 was 
detected with a TS of 17504, with a 100 MeV - 100 GeV 
photon flux of 3.68x10“^ s“^ cm“^ and Pp = 2.67. A 
lightcurve was made using bins totaling 10.5 days in GTI 
time, and ordered according to TS (a total of 88 bins). 
The progressive-binning was started from the single low¬ 
est bin, with the next-highest bin continually added as 
described above until all bins were added. At each step 
the flux (or 95% upper limit) was calculated for the five 
canonical Fermi energy bins, as for PKS 0637-752. 

3. RESULTS: TESTING THE IC/CMB MODEL 

We show in Figure [6] the radio to X-ray SEDs for the 
jet of both PKS 0637-752 (left) and 3G 273 (right). For 
PKS 0637-752 we have taken both the Chandra X-ray 
and Hubble Space Telescope (HST) infrare d and optical 
data (NIGMOS, WFPG2 and ACS) from iMehta et al.l 
(|2009ll . We have also re-derived the Spitzer infrared 
fluxes for the brightest complex of knots (wk7.8, wk8.9, 
wk9.7, wklO.7), follow ing the same methods reported in 
llJchivama et al.1 (|2005l l . We have also measured updated 
radio fluxes based on a re-analysis of archival and new 
ATCA data at 4.8, 8.4 , and 17.8 GHz (first published 
in [Godfrey et al.ll2012bll . Fo r 3C 273, data is taken from 
[Jester et al.l (|2005ll2006il and [Uchivama et ^ (|2006f ) and 
references therein. 

In both figures, we consider two scenarios: the first 
combines the photometry of the brightest/nearest knots 
to the core in order to test the I C/CMB pred i ction (black 
points and lines). In 3C 273, [Jester et al.l ([20051 ) have 
already shown that only knots A and B1 have X-ray 
indices similar to their radio indices (which is required 
for IC/CMB), so our “combined knot” scenario includes 
only these two knots. For PKS 0637-752, we use all the 
bright X-ray knots just before the turn in the jet (wk7.8, 
wk8.9, wk9.7, and wklO.6) where one might assume that 
some deceleration likely takes place. We do not include 
the only other X-ray detected knot (wk5.7) because it is 
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TABLE 2 

Results of the Fermi Data Analysis 


Source 

Band 

El 

i?2 

log Freci. 

95% Limit 

Bins 

Combined Knots* 

Single Knot^ 



(GeV) 

(GeV) 

(Hz) 

(erg cm~^) 

Added 

T^ic/Cmb % Ruled 

-Idc/CMB % Ruled 








(erg s“^ cm“^) Out 

(erg s“^ cm“^) Out 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) (9) 

(10) (11) 


0637-752 

1 

0.1 

0.3 

22.6 

1.05x10-12 

29 

9.0x10-13 

92.9 

3.6x10-13 



2 

0.3 

1 

23.1 

3.17x10-13 

32 

8.8x10-13 

99.8 

3.6x10-13 

94.5 


3 

1 

3 

23.6 

1.80x10-13 

27 

7.4x10-13 

99.98 

3.2x10-13 

98.7 


4 

3 

10 

24.1 

1.43x10-13 

50 

5.3x10-13 

99.95 

2.5x10-13 

98.6 


5 

10 

100 

24.7 

2.09x10-13 

75 

2.3x10-13 

95.9 

1.3x10-13 


3G 273 

1 

0.1 

0.3 

22.6 

2.72x10-11 

1 

2.1x10-12 


1.6x10-12 



2 

0.3 

1 

23.1 

4.63x10-12 

2 

2.8x10-12 


2.1x10-12 



3 

1 

3 

23.6 

8.20x10-13 

5 

3.6x10-12 

>99.99 

2.8x10-12 

>99.99 


4 

3 

10 

24.1 

4.46x10-13 

31 

4.5x10-12 

>99.99 

3.5x10-12 

>99.99 


5 

10 

100 

24.7 

3.56x10-13 

30 

5.2x10-12 

>99.99 

4.1x10-12 

>99.99 


Combined Knots are wk7.8, wk8.9, wk9.7, and wklO.6 for PKS 0637-752, knots A and B1 for 3C 273 
^ Single Knots are wk8.9 for PKS 0637-752 and knot A for 3C 273 


not consistently detected at other wavelengths. The sec¬ 
ond scenario assumes that the X-rays from the weaker 
knots are already not from IC/CMB, and so only the 
photometry of the X-ray brightest knot is plotted: knot 
A in 3C 273 and wk8.9 in PKS 0637-752, in both cases 
plotted as blue points and lines. The thin solid lines 
through radio-optical points show a (phenomenological) 
synchrotron spectrum fitting the data, while the heavy 
line shows the corresponding IC/CMB curve to match 
the X-ray flux levels. As shown, for both jets, the 95% 
upper limits in several bands violate the IC/CMB pre¬ 
dictions under either scenario. 

We report in Table [2] a summary of the Fermi data 
analysis for PKS 0637-752 and 3C 273. We list the defi¬ 
nition of the energy bins in columns 2-5, followed by the 
deepest 95% upper limit flux level (in vFi,) reached in 
our progressive binning for each energy bin in column 6. 
The corresponding number of bins co-added to reach that 
limit is given in column 7. In column 8 we list the flux 
predicted under the IC/CMB at the frequency given in 
column 5. This flux corresponds to the IC/CMB model 
prediction for the combination of knots wk7.8, wk8.9, 
wk9.7, and wklO.7 in PKS 0637-752, and knots A and B1 
for 3C 273. In column 10, we have calculated at what sig¬ 
nificance level our observations rule out the level of pre¬ 
dicted IC/CMB flux given in column 9. For the final two 
columns, we also give the predicted flux under IC/CMB 
and the significance-level that we can rule it out, but 
only for the X-ray brightest knot of each jet (wk8.9 and 
knot A, respectively). As shown, the IC/CMB model is 
ruled out at a > 99.9% level for PKS 0637-752 and at > 
99.99% level for 3C 273 under the combined knot scenar¬ 
ios, and ruled out at the 98.7% level for PKS 0637-752 
and at > 99.99% level for 3C 273 when considering only 
the single brightest knot. 

4. DISCUSSION 

These two cases where the IC/CMB origin for the X- 
rays has been unambiguously ruled out join with that of 
PKS 1136-135, where high UV polarization has shown 
that the second component (UV to X-ray) must be syn¬ 
chrotron in origin, since significa nt polarization is n ot ex¬ 
pected in the IC/CMB scenario (|Cara et al.ll2013ll . The 


UV polarization method is unfortunately not able to be 
applied in general, as not all of the known quasar jets 
show the second component already dominating in the 
UV, being instead dominated by the radio-optical syn¬ 
chrotron component, though we note that 3C 273 could 
be confirmed in this way. 

With IC/CMB ruled out in these three cases, we 
must explore alternative sources of X-ray emission for 
these knots. Synchrotron self-Compton (SSC) emis- 
sion was ruled out v ery early on as being far too weak 
(j Chart as et al.l[20()nl) . unless very low (far from equipar- 
tition) magnetic fields were adopted; however, the total 
power requirement then becomes far too high (> 10^® erg 
s“^) and the Doppler factor very s mall (<< 1) suggesting 
an unrealistically de-be amed jet (|Tavecchio et al.l [20001 
iDermer fc AtovanI 120041) . A viable alternative, though 
requiring further exploratio n, is a hadronic o rigin for the 
X-rays. Such models (e.s. lAharonianl l2002h would not 
be in conflict with t he high polarization measurements 
of iCara et al.l (|2013ll . and these models may be tunable 
enough to m atch the X-ray flux level while avoiding the 
GeV limits (| Zhang fc Bottche"^ 120131: iKundu fc Guotal 
l2014fl . A (leptonic) synchrotron origin for the X-rays 
from a second population of electrons is also not in con¬ 
flict with any of the data in hand, and further, relaxes 
many of the ‘uncomfortable’ constraints of the IC/CMB 
model. Very small angles to the line of sight are not 
required, and the total jet power required is consider¬ 
ably less (jPermer fc AtovanI 120041 1. as the electron en¬ 
ergy distribution need not be extended to very low val¬ 
ues. The main objection to a second synchrotron com- 
pone nt heretofore ha s simp ly been its unexplained na¬ 
ture; I Schwartz et al.l (1200011 notes that there is no rea¬ 
son why a second population of high-energy electrons 
should be co-spatial with the first. However, this co¬ 
occurrence of two very different electron populations, if 
the correct interpretation, is obviously a very important 
clue to the particle acceleration mechanism in large-scale 
jets, of which we still know little. An additional charac¬ 
teristic of the second electron population, at least for the 
case of PKS 0637-752, is that it requires a low-energy- 
cutoff in the electron energy distribution at TeV ener¬ 
gies, otherwise the spectrum should extend to optical 
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Fig. 7.— (Note: Figure description below applies to both panels; PKS 0637-752 is shown at left and 3C 273 is shown at right). The core 
flux points are shown as gray circles with phenomenological SED fit through the points drawn as a thin gray line. The big blue bump, 
visible only in 3C 373, is shown as black dashed line and not included in the beamed emission fits. In comparison, the total kpc-scale jet 
flux shown as gray triangles. The flux scale at right only applies to solid curves in the figure. The impression of the total dominance of the 
core flux is mainly a product of the beaming difference between the two components, as seen when beaming-corrected (angle-integrated) 
luminosities are plotted, rather than those assuming isotropy. Blue dashed line is the core fit times 1/5^ with (5 = 15 for both jets, 
representing the angle-integrated total core output. The gray zone is the range of possible angle-integrated output for the knots given our 
current constraints on <5 for the knots. Finally, the light blue shaded area is the range of IC/CMB emission from the knots possible given 
the same S constraints. Note that the TeV emission in particular is already constrained to be much higher than is typical for TeV blazars 
(Ri erg s“l). 


energies, contr ary to the cutoff se en in the radio optical 
knot emission (|Mehta et al.ll200^ . Acknowledging that 
further work may be warranted in the direction of evalu¬ 
ating hadronic models, we focus the rest of the paper on 
the implications for jet physics if the X-ray flux in quasar 
jets is synchrotron emission from a separate, high-energy 
electron population. 

Assuming that the X-ray emission from the jets of 
3C 273 and PKS 0637-752 is synchrotron in origin, an 
interesting consequence follows for our accounting of the 
large-scale-jet contribution to various backgrounds, espe¬ 
cially at TeV energies. Jet one-sidedness clearly indicates 
that the kpc-scale jets are at least mildly relativistic, and 
thus IC/CMB emission must occur at some level. Due 
to the very low background in the highest-energy Fermi 
bands, the flux limits reachable by Fermi's sky-scanning 
mode of operation should allow us to eventually either de¬ 
tect this emission or put very strong limits on the factor 
of B/6 which characterizes the flow on the kpc scale. The 
current Fermi 95% upper limits already constrain J 7.8 
for 3C 273 and J < 6.5 PKS 0637-752, under the assump¬ 
tion of equipartition magnetic fields, w here we take B5 
= 1.5 x 10-’^ G for PKS 0637-752 from iTavecchio et~al 
(200C) and i3i5=1.0x 10“'* for 3C 273 from iJester et al 
(j2006D . Lower magnetic field values would only decrease 
these upper limits on 5. These limits are already low 
enough to have interesting consequences for our under¬ 
standing of the total radiative output of AGN jets on the 
kpc scale. 

It is generally assumed that the radiative output of 
quasar jets is dominated by that occurring at the ‘core’, 
the base of the jet which is presumed to be very near the 
black hole (or ^ parsecs away at most) and is therefore 
unresolved even in VLBI imaging. Certainly, the ob¬ 
served fluxes are dominated by this part of the jet due to 
Doppler beaming whenever the jet is pointed fairly along 
our line-of-sight. This is depicted in Figure [71 where the 
core flux of both jets are shown as gray circles, and the 


total kpc-scale jet as gray triangles (flux scale on right 
axis). The luminosity scale at left applies to these points 
only under the incorrect assumption of isotropy. If we 
correct these values for beaming, we get the real ‘angle- 
integrated’ total power output from the core and the jet, 
plotted as a dark blue dashed line and the gray shaded 
area, respectively. VLBI observations of superluminal 
motions place a lower limit on r=15 for the cores of both 
3C 273 and PKS (3637-7 52, respectively ([Lister et al.l 
120131 : lEdwards et al.l 120061) . We have applied a correc- 
tion (multiplying by l/J^. iGhisellini fc Tavecchioll20inD 
assuming J=15 for the cores of both jets to give 
the angle-integrated core luminosity (dark blue dashed 
line). To calculate the angle-integrated luminosity of the 
knots, we apply a lower limit value of J = 1.9 which 
comes from statistical argu ments based on populations 
([Arshakian fc Longaiill200J ). which gives the dotted-line 
upper edge to the gray shaded area, while the current 
5 limits from Fermi give the lower dashed-line limit. 
The true angle-integrated luminosity of the knots is thus 
somewhere in the gray zone. 

Gomparing the dashed-blue curve and the gray shaded 
area, it is interesting to note that the knots are appar¬ 
ently not necessarily insignificant in total output when 
compared to the core. The jet power dissipated by ra¬ 
diative losses appears to occur over much larger physical 
scales than has previously been appreciated. Definite 
conclusions will require tighter constraints on the 5 fac¬ 
tors of both the core and the knots, but it is possible that 
large-scale jets contribute more than the core in the UV 
to X-rays, in addition to their general dominance in the 
radio. Large-scale jets could thus also be a contributor 
to some astrophysical backgrounds. 

A further observation follows from the realization that 
the X-rays are synchrotron in origin: the electrons pro¬ 
ducing the synchrotron X-rays will themselves upscat- 
ter the GMB to produce a GeV to TeV spectrum. The 
angle-integrated total power in the IC/CMB component 
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is shown in Figure [7] as a light blue shaded area. Note 
that the bounds in this case are flipped; the upper d 
limit forms the upper edge of the allowed zone. This 
is because IC/CMB emission increases with the Doppler 
factor as 6^ even after the angle-integrated 47r luminosity 
correction is included. Even in the minimum i5=1.9 case, 
these jets are already constrained to produce fluxes in ex¬ 
cess of 10^^ erg s“^ which is the typical angle-integrated 
luminosity for canonical TeV blazars (adopting an ob¬ 
served luminosity of erg s“^ and a modest (5=10). 

This is important because it has been proposed that 
(lower-power) TeV BL Lacs are signihcant contributors 
to the heating of the intergalactic medium (IGM; e.g. 
Bro(lerick et ai.l[2012l:lChang et al.ll2012l:lLamberts et al.l 

20151 but see also iSironi &: GianniosI l2014f) . where the 
high-energy gamma-rays (> 100 GeV) produced by these 
jets pair-produce off the extragalactic background light, 
eventually depositing kinetic energy through plasma in¬ 
stabilities. Clearly, if this mechanism is important, the 
large-scale jets of powerful quasars may be a far more im¬ 
portant source of > 100 GeV photons than TeV BL Lacs, 
and possibly the dominant class of sources for blazar 
heating. 

We have not applied an EBL correction to these spec¬ 
tra merely to illustrate the total intrinsic output. How¬ 
ever, EBL absorption is very important at TeV energies, 
and would make direct observation of this TeV compo¬ 
nent difficult. Even assuming the most optimistic case of 
^=7.8 for 3C 273. at redshift 0.158 the EBL absorption 
(jFinke et al.ll201fll l is already high enough that it would 
take at least 100 hours of observations by the future CTA 
mission to detect the beamed IC/CMB component at 
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TeV energies. Thus it is unlikely that many anomalous 
X-ray jets will have a synchrotron origin for the X-rays 
directly confirmed via TeV observations, though if Fermi 
begins detecting the IC/CMB component, an upturn at 
the highest energies might be visible in a few cases. The 
remaining best direct observation is via polarization - 
either in the UV, for those that show the second compo¬ 
nent emerging there, or with future X-ray polarimeters. 
Finally, we note that as long as Fermi continues to op¬ 
erate, the low background at the highest energies should 
allow continually improving constraints on the 5 factors 
of large-scale jets. 

5. CONCLUSIONS 

We have shown that the expected gamma-ray emission 
required if quasar jets produce X-rays via IC/CMB has 
been ruled out in two quasar jets, 3C 273 and PKS 0637- 
752, at a high significance level using Fermi upper limits. 
Examining the evidence, we favor a synchrotron origin 
for the X-rays of their large-scale jets. This considerably 
relaxes the power requirements for the jet away from 
the near and super-Eddington values. We have shown 
that the limits on the IC/CMB gamma-rays constrain 
the jets to have 5 values of <7.8 for 3C 273 and <6.5 for 
PKS 0637-752 assuming equipartition. A very interest¬ 
ing outcome of the synchrotron explanation of the X-rays 
is the prediction that quasar jets may radiate far more at 
TeV energies than TeV BL Lacs. Future work remains 
to see if 3C 273 and PKS 0637-752 are outliers, or if 
IC/CMB can be ruled out as the source of the anoma¬ 
lous X-rays for more quasar jets. 


Kharb et al. 2012, ApJ, 748, 81 

Kundu, &; Gupta. 2014, MNRAS, 444, L16 

Lamberts et al. 2015, ArXiv e-prints 

Lister et al. 2013, AJ, 146, 120 

—. 2009, AJ, 138, 1874 

Marshall et al. 2011, ApJS, 193, 15 

—. 2005, ApJS, 156, 13 

Mehta et al. 2009, ApJ, 690, 1706 

Meyer, & Georganopoulos. 2014, ApJ, 780, L27 

Mullin, & Hardcastle. 2009, MNRAS, 398, 1989 

Nolan et al. 2012, ApJS, 199, 31 

Sambruna et al. 2004, ApJ, 608, 698 

—. 2002, ApJ, 571, 206 

—. 2001, ApJ, 549, L161 

Schwartz et al. 2000, ApJ, 540, L69 

Siemiginowska et al. 2003, ApJ, 598, L15 

—. 2007, ApJ, 657, 145 

Sironi, &: Giannios. 2014, ApJ, 787, 49 

Tavecchio et al. 2000, ApJ, 544, L23 

Uchiyama et al. 2006, ApJ, 648, 910 

—. 2005, ApJ, 631, L113 

Wilson, & Yang. 2002, ApJ, 568, 133 

Zhang, & Bottcher. 2013, ApJ, 774, 18 


ETM acknowledges Fermi Grant NNXI3A088G. 

Facilities: Fermi, HST (ACS, NICMOS, WFPC2), CXO (ASIS), Spitzer, ATCA. 











